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ABSTRACT 

We present threshold enhanced QCD corrections to the bottom quark energy distribution in Higgs 
boson decay and to hadroproduction in Z+Z annihilation beyond leading order in the strong cou- 
pling constant. This is achieved using the resummed decay distribution obtained using renormal- 
isation group invariance and the mass factorisation theorem that they satisfy and Sudakov resum- 
mation of soft gluons. 



The Higgs boson (H) which emerges in the electroweak symmetry breaking mechanism [1] is 
the only particle that is yet to be discovered in the Standard Model(SM). Searches at LEP exper- 
iments [2] indicate that the lower bound on its mass is around 114.2 GeV and the upper bound 
around 260 GeV at 95% confidence level. The hadronic machine Tevatron at Fermi-lab currently 

running with increased energy and the upcoming Large Hadron Collider (LHC) at CERN have the 
physics goal of discovering the Higgs boson. At hadronic colliders, among the potential channels 
through which the Higgs boson can be produced, the vector boson fusion process is one of the 
promising ones [3]. Here the Higgs boson is produced in a reaction Pi+ P2 ^ V + H, where Pi 
are the incoming hadrons, V is the vector boson W^/Z, followed by the dominant decay mode 
H ^b + b. For the heavy Higgs boson searches, the decay to bottom quarks can be accessible in 
the Two-Higgs Doublet models due to enhanced couplings [4]. 

A similar process in nature is inclusive hadroproduction in annihilation [5]. Here the 
incoming leptons annihilate into vector boson such as y* /Z which then fragments into hadrons, 
i.e., — > Y*,Z — > P+X, where P is any hadron and X is the remaining final state. The double 
differential cross section for producing a hadron with energy fraction x of the parton produced in 
the annihilation is given by 

= -(1+cos^e) — ^(jc,^^) + -sm^e — ^(;c,^^) + -cose — ^(jc,^^) , (1) 



JxJcosG 8 dx A dx A dx 

where is the center of mass angle of the final state hadron. The energy fraction is given by 
x = lp-q/q^, where p and q are the momenta of the final state hadron and the intermediate vector 
boson respectively. The B hadron production in channel is an important process to determine 
the non-perturbative b quark fragmentation functions [6]. Precise knowledge of them will reduce 
theory uncertainties in Higgs decays to bottom quarks and also in the top quark mass reconstruc- 
tions in the top quark decays at Tevatron and LHC. 

In the case of the decay process H ^ b + b, we study the energy distribution of one of the 
bottom quarks. We use the perturbative fragmentation approach [7] which is valid when nif, <^ nin 
to factorise the decay distribution into a part which contains the Higgs boson decay to massless 
partons and a part containing perturbative fragmentation functions denoted by that describe 
the fragmentation of massless partons into massive bottom quarks. In this approach, the large 
logarithms log(m^/m|), where mn is the Higgs boson mass and bottom quark mass, can be 
resummed using perturbative fragmentation functions. The normalised decay distribution is given 
by 



1 dFb 



{xb,q^,ml) = [ —Ciiz^l^^l^^l^F)^! (—,ml,iJ^^ I = b,g, (2) 

J Jxb Z \ Z / 



To dxh 

where Tq is the Born decay distribution, X}, = Ipi, ■ q / with p/,, q the momenta of b quark and the 
Higgs boson, with masses m^, and mn respectively. Here q^ = mj^. The renormalisation scale jjr is 
due to ultraviolet renormalisation and the factorisation scale /jp is due to mass factorisation of the 
coUinear singularities. The fixed order result for Cj for I = b up to next to leading order (NLO) 
in strong coupling constant is available in the literature [8]. The soft gluons that result from 
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the outgoing partons lead to large logarithms that can be resummed systematically. The impact of 
such effects up to the next to leading logarithmic (NLL) level has been studied along with fixed 
order NLO QCD corrections in [6, 8]. 

The inclusive hadroproduction in /+/ annihilation process factorises as 



where z = Ipi -qjc^ \s, the partonic energy fraction. Here, pi is the momentum of the parton that 
fragments into the final state hadron with momentum p. The non-perturbative function (x/z, 
is the fragmentation function that describes the fragmentation of the parton of type / into hadron 
of type P. The perturbatively calculable coefficient function C/ is known up to NNLO in QCD 
for 7 = q,q,g, [9-12]. The NNLL soft gluon resummation to this process is also known in the 
literature [13-17]. 

Even though perturbative Quantum Chromodynamics (pQCD) provides a framework to suc- 
cessfully compute various observables involving hadrons at high energies with less theoretical 
uncertainty for the physics studies, the fixed order QCD predictions often have limitations in ap- 
plicability due to the presence of various logarithms that become large in some kinematic regions 
which otherwise can be probed by the experiments. The standard approach to probe these regions 
is to resum the class of such large logarithms supplemented with fixed order results. This can 
almost cover the dominant kinematic region of the phase space. 

In this paper, we follow the method that we used in Refs. [18, 19] to study the soft gluon effects. 
This generalises our earlier approach to include any infrared safe decay distribution in perturbative 
QCD. We systematically formulate a framework to resum the dominant soft gluon contributions 
in z space to the decay distributions, where z is the appropriate scaling variable that enters in the 
process. We have followed a similar approach described in [18, 19] that uses renormalisation group 
(RG) invariance, mass factorisation and Sudakov resummation of soft gluons as the guiding princi- 
ples. Using the resummed results in z space, we compute the soft plus virtual part of the dominant 
partonic decay distributions beyond NLO. In Ref. [18] we determined the threshold exponents 
and Bi ois up to the three loop level for Drell-Yan process, Higgs boson productions and deeply 
inelastic scattering (DIS) cross sections using our resummed soft as well as soft plus jet distribu- 
tion functions. Here we extend the similar all order proof which establishes the relation between 
soft plus jet distribution functions and the threshold resummation exponents and demonstrate the 
usefulness of this approach to derive higher order threshold enhanced corrections for any infrared 
safe decay distributions. 

We find that the soft plus jet distribution relevant for our present study can be gotten entirely 
from that of DIS due to the crossing symmetry between them [20-24]. In fact, we find that they are 
identical. In addition, the Altarelli-Parisi(AP) splitting functions that determine the scale evolution 
of the fragmentation functions coincide with the AP splitting functions of parton distribution func- 
tions in the threshold region. Extensive discussions on this topic can be found in [25-30]. With 
these two essential ingredients and using the method described in [18, 19] we could successfully 



1 d<3j{x,q^) 




i = q,q,g 



(3) 



a(0) dx 
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reproduce the fixed order NLO soft plus virtual part as well as the NLL resummation exponents [8] 
for the bottom quark production in Higgs decay and also the fixed order NNLO soft plus virtual 
part [9-12] and the NNLL resummation exponents [13-17] for the hadron production in an- 
nihilation. In addition, we can compute the fixed order NNLO soft plus virtual part of the b quark 
energy distribution in the Higgs boson decay and also partial but dominant N^LO results for / = 3, 4. 
With the same level of accuracy, we can obtain N^LO with i = 3,4 soft plus virtual contributions 
to hadroproduction in /+/~ annihilation. 

In this article, we will mainly concentrate on dominant contributions coming from the thresh- 
old effects. This contribution to the energy distribution of the bottom quark in Higgs boson decay 
and to the hadroproduction in the Z+Z^ annihilation process can be obtained by adding the soft 
part of the decay distributions with the renormalised virtual corrections and performing mass fac- 
torisation using appropriate mass factorisation counter terms. We call this infra-red safe combina- 
tion a "soft plus virtual"(5v) part of the decay distribution. The soft plus virtual part denoted by 
A*)'(z, q^iiJ^ii^) of the perturbative decay distributions C/ after mass factorisation is found to be 



(4) 

8=0 



where ^^(z,<7^,/j|,/j^,£) is a finite distribution. For a Higgs boson decaying to bottom quarks, 
we put I = h and for that decaying to gluon, I = g. For the hadroproduction in Z+Z^ annihilation, 
/ = q,q. Here {z.,q^,iJi\,i^,2) is computed in 4 + 8 dimensions, where n denotes the number of 
space-time dimensions. 

^'\z,q^AA.^) = (^ln(z^(4,;"|,;"',e))' + ln|F^(a„(22,/|2 8)|2^5(l-z) 

+2$^(4,^^^2 z,e)-Clnf//(4,^^^^,z,E), I = b,g,q,q (5) 

The symbol "C" means convolution. For example, C acting on an exponential of a function f{z) 
has the following expansion: 

= 5( 1 - z) + + ® /(z) + i/(z) ® /(^) ® /(^) + . . . (6) 

The function f{z) is a distribution of the kind 5(1 — z) and 'Dj, 

"In'(l-z) 



/ = 0,1,--- (7) 



and the symbol denotes the Mellin convolution. We drop all the regular functions that result 
from various convolutions. F\ds,Q^,i-^,z) are the form factors that enter these processes. For 
the Higgs boson decay to bottom quarks, it is related to the Yukawa type of interaction and to the 
gluons, it is related to Higgs boson field coupled to the kinetic term of the gauge fields. For the /+/~ 
annihilation, the form factor arises from the vector current. In the form factors, = —q^ = — 
for the Higgs boson decay and = —q^ = —q'?^/z- functions ^ \ds.,q^,iiP',z,&) are called 
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the soft plus jet distribution functions. The unrenormalised (bare) strong coupling constant ds 
is defined as — gl/{l6'!i^) where gs is the strong coupling constant which is dimensionless 
in n = 4 + £. The scale ^ comes from the dimensional regularisation in order to make the bare 
coupling constant gs dimensionless in n dimensions. The bare coupling constant ds is related to 
the renormalised one by the following relation: 



Sea, = Z{nn)as{^R) y-^j (8) 

The scale is the renormalisation scale at which the renormalised strong coupling constant 
as{ij\) is defined. The factorisation scale fjp is due to mass factorisation and the constant 5e = 
ejcp{(y£ — ln47r)£/2} is the spherical factor characteristic of n-dimensional regularisation. The 
renormalisation constant up to three loop is given by 

= l+as{pR) — + a,{/UK)\ + — +a.W) -.0 +^ (9) 



£2 e y ' -syrKjy g3 3^2 3e ^ 

The renormalisation constant Z(/j|) relates the bare coupling constant ds to the renormalised one 
as{fj\) through Eq. The coefficients Po^Pi and P2 are the coefficients that appear in RG equa- 
tion for the strong coupling constant up to three loop level [31, 32] 

The factors Z^(ai,/j|,/j^,£) are the overall operator renormalisation constants which renor- 
malise the bare form factors F\ds, Q^,/j^,e). For the vector current Z^{ds.,id\.,i^ ,£) = 1, but for 
the bottom quarks and gluons coupled to the Higgs boson, the corresponding form factors get 
overall renormalisations [31, 33, 34]. The bare form factors F\ds, £) corresponding to the 

unrenormalised operators satisfy the Sudakov differential equation [35-38]. In dimensional reg- 
ularisation, the formal solution to the differential equation up to four loop level can be found 
in [18, 39] in terms of A\ the standard cusp anomalous dimensions, and the constants G/(£) for 
both I = q,b and I = g [40] to the required accuracy in £. The single poles of the form factors 
contain the combination [18, 19,39,41] 



2(^5/-5/,,/|3,_i-5/,,7f„i)+// 



at order aj,. The terms proportional — 2(5/^g / |3/_i +5/i,7f_j) come from the large momentum 
region of the loop integrals that are giving ultraviolet divergences. The poles containing them 
will go away when the form factors undergo overall operator UV renormalisation through the 
renormalisation constants Z ' which satisfy the RG equations 

d 



J CO 

^l-^\nZ\ds,^W,t) = l^aM)y^_i (10) 

where £ ^ is set. The constants / |3;-i and are anomalous dimensions of the renormalised 
form factors and F^ respectively. After the overall operator renormalisation through Z ^ and 
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coupling constant renormalisation through Z, the remaining poles will be proportional to and 
fi^ in addition to the standard cusp anomalous dimensions A/. The constants B\ and f- are known 
up to order [40^3]. 

The coUinear singularities resulting from massless partons are removed in MS scheme using the 
mass factorisation kernel T{z,i^,2). The kernel r(z,jU^,e) satisfies the following renormalisation 
group equation: 

/i|.-^f(z,/i2,e) = ip(z,//2)^f (z,/i2,e) . (11) 
The function P{z,i^) are the well known (matrix- valued) Altarelli-Parisi splitting functions 

oo 

nzA) = 'L<{i/p)p^'-'\z) (12) 

The diagonal terms of the splitting functions P'^^ (z) have the following structure 



5|+i5(l-z)+A|+i©o 



where P^^^ jj are regular when z — > 1. In the case of the soft plus virtual part of the decay distribu- 
tions, only the diagonal parts of the kernels contribute. In the MS scheme, the kernel contains only 
poles in e. The kernel can be expanded in powers of bare coupling ds as 

f(z,/|2,e) = 5(l-z) + f;4f4) '4f«(z,e) • (14) 

i=\ \j" / 

The constants r(') (z, e) are expanded in negative powers of e up to the four loop level can be read 
from Ref. [18] where similar RG equations were solved. 

It is natural to expect that the soft plus jet distribution functions have a pole structure in 8 
similar to that of and F// so that the decay distributions A^J are finite in the limit £ — 0. This 
implies that they satisfy a Sudakov type differential equation that the form factors do. Solving 
the Sudakov differential equation for ^\ds,q^,i^, z, e), we get 

$^(fl„^^//2^Z,8) = ^\a,,q^{\-z),l?,^) 



!=1 

where $ coincide with the that appear in the deep-inelastic scattering cross section 

[19]. This is the result of exact crossing symmetry between the processes under study and the DIS 
process in the threshold region. The constants ^^j^^ (£) ior I — q corresponding to DIS are given 
in [18, 19]. Since the b quark is treated massless, we use the same constants for the Higgs boson 
decay also. 
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The threshold corrections that dominate when the partonic scaling variable z approaches its 
kinematic limit, which is 1, through the distributions 6(1 — z) and can easily be resummed 
in Mellin because the decay distributions appear as convolution of partonic distributions and 
fragmentation functions. This has been a successful approach due to several important works 
(see [13-17]). The higher order threshold exponents can be found in [18,44-46]. We find that 
the soft plus jet distribution function {ds,q^,iJ^,z,£) captures all the features of the N space 
resummation approach by expressing (fTSt as 



oo / 2 \ 't 

+5(1 -z) S'^^^f\e) 



,2X.1 



+ (.,0^)^ E^;(^) ^e^'"'W C"^) 



where 



Gl{aM\l-z)).z) = ia'J'-^^\'s^^G^f\e). (17) 

i=l V A* / 

Here, the constants Gsf'\e),K ^'^'\e) and ^sf\e) can be found in [18, 19]. The second line of 
the equation ([T6b contains the right poles in e to cancel those coming from the form factor as well 
as from the mass factorisation kernel. It also contains the terms that are finite as £ becomes zero 
through the constants g jg ■(£) that appear in G^f'\£) (see [19]). Since they are multiplied by 5(1 — 
z), the fixed order soft plus virtual part of the decay distributions gets contributions from them at 
higher orders. The third line in the eqn. (fT6t contains only poles in £ and they all cancel against ©o 
parts of the mass factorisation kernel. Hence, adding the eqn. (fT6b with the renormalised form factor 
and the mass factorisation kernel and performing the coupling constant renormalisation and then 
finally taking Mellin moment, we reproduce the exponents in resummation formula(see [13-17]) 
(after setting £ ^ 0). We find that the function G^j (a^ (^^(1 — z)) ,£) appearing in the first line of 
eqn.(fT6t coincide with the exponent 5^^^.^^ {as [q^{l—z))) appearing in the resummation formula 

oo 

BLay{^s{qHi-z))) = I, <. {q\\ - z)) B^,^^^, 



!=1 



Gi(a, {qHl-z)).£) 



(18) 

£=0 



In addition, we get a new exponent that comes from the Mellin moment of 5( 1 — z) part and hence 
A'^ independent. The Mellin space resummation exponents contain the cusp anomalous dimensions 
and the constants B'^^^^^ ■ which we find are identical to those of DIS. These constants are known 
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up to three loop level whose numerical impacts on the physical observables will be important in 
order to reduce theoretical uncertainties coming from renormalisation and factorisation scales. 

Using the z space resummed expression given in Eq. Q and the known exponents, we present 
here the results for A"'*^'^ for the bottom quark energy distribution in Higgs boson decay beyond 
leading order in QCD. We present the complete soft plus virtual contribution to order N^LO for 
/ = 1,2 and a partial N^LO result, i.e., without 5(1 —z) part. In addition, we extend this approach 
to N^LO order where we can obtain the partial soft plus virtual contribution coming from all 
except J = 0, 1 for the N^LO coefficient A^^'*^^\ Here also we cannot determine the 5(1 —z) part. 
We repeat this for the hadroproduction in /+/ annihilation up to N^LO order. These fixed order 
results are expected to reduce the renormalisation and factorisation scale uncertainties. Expanding 
AYiz,q^,lul,fjf) in powers of as{ijl) as 



i=0 



and choosing ij^= /Jp = for simplicity, we find 



r.v,(0) 

:.^v,(i) 



:sv,{2) 



= 5(1 -z) 
= 5(1 -z) 

= 5(1 -z) 



Cf 3 + 8^2 



+ ©0 



Cf 



+ ©1 



Cf 4 



(19) 

(20) 
(21) 



HfCF -91/12-14/3 C2 + 4/3 ^3)+CfCa 1691/24 + 95/3 ^2 



-49/5 C2 + 32/3 ^3)+Cj (109/8 + 31 ^2 + 30 ^^-78 ^3 
+©0 rifCF (247/27- 8/3 C2)+CfCa (-3155/54 + 44/3^2+40^3) 
+C^ (-21/2-8^3) +©1 h/Cf (-58/9) +CfCa (367/9-8^2) 
+C^ (21 + 16^2) +2)2 n/Cf (4/3) +CfCa (-22/3) +C| (-18 



+©3 



Cj (8 



(22) 



©0 



HfCFCAll 60906/729 - 9920/8 1 ^2 + 208/ 1 5 - 776/9 ^3 



+nfC} (l6423/108- 722/27 ^2 - 16 C2 - 60 ^3) +«^ (-8714/729 
+232/27 ^2-32/27 ^3) +Cf Cj (-599375/729- 176/3 ^2 ^3 + 32126/81 C2 
-652/15 ^^ + 21032/27 ^3 -232 ^5) +C^Ca (-31151/72 + 80^2 ^3 
+3365/27 C2 + 39 + 1988/9 ^3 - 120 ^5) + C| ( - 479/8 - 64 ^2 ^3 

+ ©1 h/CfCa (- 15062/81+512/9 ^2 



-45 ^2 + 6^2 + 178 ^3+432 ^5 
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+16 ^g) +n/C| (I - 1343/9 + 64/3 ^2+ 112/3 Cs) +«/ Cf (940/81-32/9 ^2) 
+Cf Cl (50689/81-680/3 ^2+ 176/5 ^2-264 ^3) +Cj Ca (i3783/18 

-352/3 C2 - 196/5 -592/3 ^3) +Cj (l81/2 + 4 ^2 - 104/5 - 360 ^3) 

HfCpCA (1552/27 -16/3 C2)+n/C| (827/9-64/3 ^2) 
+n}CF (-II6/27) +CfC] (-4649/27 + 88/3 C2)+C|Ca (-10009/18 
+460/3 C2 + 240 C3) + 4 ( - 153/2 + 72 ^2 + 16 C3) 



+©2 



+©3 



HfCpCA (-176/27) +n/C| (-28O/9) +n}CF (16/27) 



+Cf Cl (484/27) +Cj Ca (i732/9 - 32 C2) +4 (60) + ©4 4 (40/9 
+C| Ca ( - 220/9) + 4 ( - 30) 



+ ®5 



4 (8 



(23) 



2)2 



HfCpCl (17189/9-5096/9^2+176/5 C2- 352^3) 

+n/ Ca (964334/243 - 57524/27 ^2 + 2332/15 C,l - 1 1032/9 ^3) 
+nf Cl (12299/9 - 2336/3 ^2 - 728/5 C2 - 936 ^3) +«/ Cf Ca ( - 7403/27 
+688/9 ^2 + 16 ^3) + n) Cj 7 1776/243 + 4088/27 ^2 + 304/9 ^3) 
+n}CF (940/81-32/9C2)+CfC1 (-649589/162 + 4012/3 C2 -968/5 ^2 
+1452 ^3) + C| C] ( - 6034493 /486 - 832 ^2 ^3 + 63764/9 ^2 - 2450/3 
+25336/3 C3 - 1392 ^5) +c|^ Ca ( - 12529/3 - 768 ^2 ^3 + 13718/3 ^2 
+1094/5 C2 + 13300/3 ^3 - 720 Cs) + C^ ( - 420 - 448 ^2 ^3 
+ 198 C2 + 324 + 1584 ^3 +4128 Cs) 



+2)3 



n/CFCl (-9502/27 + 352/9^2) +n/C|CA (-358142/243 

+11984/27 C2 + 1216/9 ^3) +nfCj (-5660/9 + 2608/9 ^2+ 1888/9 ^3' 
+n} Cf Ca (1540/27 - 32/9 ^2) +«/ C| (25966/243 - 736/27 ^2) 
+n} Cf (-232/81) +Cf C] (55627/81-968/9 C2) +Cj Cl (2259107/486 
-47104/27 C2 + 864/5 Cl - 13024/9 ^3) +C|- Ca (26720/9 - 16192/9 C2 
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+248/5 Ci- 11392/9 ^3) (^533/2-376 C2- 1488/5 C2- 1072 ^3) 

+©4 «/ Cf Cl (^242/9) +n/ C| Ca (^8120/27 - 80/3 ^2) +«/ C|- (^6070/27 

-560/9 C2) +«/ Cf Q (-44/9^ +4 C} (-640/27^ (8/27^ 

+Cf C] - 133 1 /27) + cl - 24040/27 + 440/3 ^2) + C| Ca ( - 35755 /27 

+4160/9 ^2 + 400 ^3) + 4 ( - 150 + 240 C2 + 400/3 Cs) 



+2)5 



n/ C| Ca ( - 704/27) + n/ C^ ( - 164/3) + C| (^64/27 
+C| cl (1936/27) +C^ Ca (998/3-48 C2) +C^ ^78 - 32 ^2) 



+©6 



UfCl (56/9) + C| Ca ( - 3O8/9) + C^ ( - 28 



+ !»7 



C 



I6/3) 



Similarly for hadroproduction in /+/ annihilation, we find 



T*v,(0) 


T*v,(0) 
^9 


= 




A^v,(l) 


X^v,(l) 


= 8(1 


-z) 


'^b 




= S(l 


-z) 






+c 


K- 


^b 




= !»0 


UfC, 



Cf(i2) 

n/CF ( - 365/ 1 8 + 8 C2) + CfCa (5269 /36 - 40 C2 - 36 ^3) 



(24) 

(25) 
(26) 



+ ©0 



C|( -36 



+ ®i 



^2 
Cf 



(4S) 



(27) 



n/C|(3071/18-56 C2) +C|-Ca(- 41047/36 + 296 C2 + 588 C3) 



+ 2)1 



n/C|(- 1426/9 + 32 ^2) 



'r^v,(4) riv,(4) 



+©2 
©2 



+ 2>3 



Cf 



+C^ (261/4 + 78 C2-96 C3) 

+C|Ca (9673/9 - 256 ^2 - 144 ^3) + 4 ( - 3 + 88 C2) 

n/C^ (16) + cICa ( - 88) + C^ ( - 216) 

UfClCA (279O8/27 - 176 ^2 -48 C3) +«/C|- (4148/3 - 1304/3 ^2) 
+n}c|(- 2122/27 + 32/3 ^2) +c|cl (- 169535/54+ 1936/3 ^2 
+264 C3) + cICa ( - 265 19/3 + 8252/3 C2 + 3528 ^3) + 4 (459/2 

n /C|Ca ( - 704/9) + n/C| ( - 4820/9 



(28) 



+1332 C2+I92 C3 



+ ©3 



+64 C2) + (64/9) + C|C1 (1936/9) + C^Ca (31322/9 



-704 ^2 -288 ^3) +4 (210 -208 ^2) +©4 
+4Ca(-880/3)+4(-360) +©5 4 (96 



n/Q 160/3 



(29) 



where the colour factors for the SU (N) gauge group are Ca = N,Cf = {N^ — 1)/ (2A^) and «/ is 
the number of active flavours. The above expressions are easily transformed to Mellin space. The 
functions "Dk are represented by polynomials of single harmonic sums only, cf. [47]. The soft re- 
summation terms supplement the corresponding representations for the 2-loop Wilson coefficients 
of the production cross sections in Mellin space derived in [48]. 

To summarise, we have systematically studied the soft plus virtual contributions to the bot- 
tom quark energy distribution in Higgs boson decay and the hadroproductions in /+/ annihilation 
processes using the formalism of threshold resummation to infrared safe decay distributions. This 
was achieved using renormalisation group invariance and Sudakov resummation of soft gluons and 
the factorisation property of these decay distributions. We have also shown how these resummed 
distributions are related to resummation exponents that appear in Mellin space. Using this ap- 
proach we have computed the soft plus virtual decay distributions at NNLO and partial results at 
N^LO and N^LO in perturbative QCD for the bottom quark energy distribution in Higgs decay and 
hadroproductions in /+/ annihilation. 
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